NASA TECHNICAL NOTE 


>© 

CO 

to 


CASE FILE 
COPY 


ANALYSIS OF LATERAL-DIRECTIONAL 
STABILITY CHARACTERISTICS OF A 
TWIN-JET FIGHTER AIRPLANE 
AT HIGH ANGLES OF ATTACK 

by Joseph R. Chambers and Ernie L. Anglin 

Langley Research Center 
Langley Station , Hampton, V a. 



NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D. C. 


AUGUST 1969 


1. Report No. 2. Government Accession No. 

NASA TN D-5361 

3. Recipient's Catalog No. 

4. Title ond Subtitle 

ANALYSIS OF LATERAL-DIRECTIONAL STABILITY CHARACTERISTICS 
OF A TWIN-JET FIGHTER AIRPLANE AT HIGH ANGLES OF ATTACK 

5. Report Date 

August 1969 

6. Performing Organization Code 

7. Author(s) 

Joseph R. Chambers and Ernie L. Anglin 

8. Performing Organization Report No. 

L-6723 

9. Performing Organization Name and Address 

NASA Langley Research Center 
Langley Station 
Hampton, Va. 23365 

1 0 . Work Uni t No. 

126-62-01-01-23 

11. Contract or Grant No. 

13. Type of Report and Period Covered 

Technical Note 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 

14. Sponsoring Agency Code 


15. Supplementary Notes 


16. Abstract 


An investigation was conducted to determine the factors producing a directional divergence at high 
angles of attack for a twin-jet swept-wing fighter airplane. The study consisted of static wind-tunnel tests, 
tuft-flow visualization tests, and calculations of the dynamic lateral-directional stability characteristics of 
the airplane. Several modifications to the basic configuration were evaluated in an attempt to delay or 
eliminate the instability. 


17. Key Words Suggested by Author(s) 

Lateral-directional stability 

Dynamic stability 

Stability at high angles of attack 


18. Distribution Statement 

Unclassified - Unlimited 


19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21. No. of P ages 

22. Price* 

Unclassified 

Unclassified 

47 

$3.00 


For sale by the Clearinghouse for Federal Scientific and Technical Information 
Springfield, Virginia 22151 







ANALYSIS OF LATERAL-DIRECTIONAL STABILITY CHARACTERISTICS 
OF A TWIN-JET FIGHTER AIRPLANE AT HIGH ANGLES OF ATTACK 

By Joseph R. Chambers and Ernie L. Anglin 
Langley Research Center 

SUMMARY 

An investigation was conducted to determine the factors producing a directional 
divergence at high angles of attack for a twin- jet swept-wing fighter airplane. The study 
consisted of static wind-tunnel tests, tuft-flow visualization tests, and calculations of the 
dynamic lateral-directional stability characteristics of the airplane. Several modifica- 
tions to the basic configuration were evaluated in an attempt to delay or eliminate the 
instability. 

The results of the investigation indicated that the directional divergence exhibited 
by the airplane was brought about by a simultaneous loss of directional stability and effec- 
tive dihedral at high angles of attack. The loss of directional stability resulted from a 
combination of an adverse sidewash region at the rear of the airplane and a reduced 
dynamic pressure at the vertical tail location. The adverse sidewash was generated by 
the wing-fuselage combination and was related to stalling of the leading-wing panel during 
a sideslip at high angles of attack. The loss of effective dihedral was also attributed to 
leading -wing-panel stall. The apparent directional divergence was determined to be, in 
reality, a highly unstable Dutch roll oscillation. The only geometric modification studied 
that significantly delayed the divergence was wing leading-edge droop. 

INTRODUCTION 

The National Aeronautics and Space Administration is currently conducting a series 
of investigations of the poststall characteristics of a high-performance swept-wing fighter 
airplane. Recently, concern has arisen over the existence of directional divergence 
(sometimes termed ’’nose slice”) at angles of attack near the stall. Tactical training and 
air combat maneuver requirements imposed on the airplane have resulted in operational 
angles of attack near the stall and the associated directional divergence which in turn has 
produced inadvertent poststall gyrations and spins. Inasmuch as instabilities of this type 
can seriously limit the maneuvering capability of an airplane, the present investigation 
was conducted (1) to identify the various factors producing the directional divergence and 
(2) to define geometric modifications or fixes which might eliminate or postpone the 


instability to angles of attack farther removed from the operational flight envelope. The 
study consisted of static wind-tunnel force tests, flow visualization tests, and calculations 
of the dynamic lateral-directional stability of the airplane. 

SYMBOLS 

All aerodynamic data with the exception of lift and drag are presented with respect 
to a body system of axes. Moment data are presented with respect to a center-of-gravity 
position of 33 percent of the wing mean aerodynamic chord. Dimensional values herein 
are given in both U.S. Customary Units and in the International System of Units. 

A,B,C,D,E coefficients of lateral-directional characteristic equation (see appendix A) 

b wing span, ft (m) 

c mean aerodynamic chord, ft (m) 

c t mean aerodynamic chord of horizontal tail, ft (m) 

Cj) drag coefficient, Fj-^q^S 

Cl lift coefficient, F^/q^S 

Cj rolling- moment coefficient, M^/q^Sb 

C m pitching- moment coefficient, My/q^Sc 

C n yawing-moment coefficient, Mz/q^Sb 

Cy side-force coefficient, Fy/q^S 

Dfo differential operator, d^ds^ 

F d drag force, lb (N) 

Fl lift force, lb (N) 

Fy side force, lb (N) 

it horizontal tail deflection (positive when trailing edge is down), deg 
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Ix 

IY 

J Z 

x xz 

k X 0 

k Z 0 

K x 0 

K Zo 

ly 

m 

M 

M X 

My 

Mz 

P 

P 

% 


moment of inertia about longitudinal body axis, slug-ft^ (kg- m2) 
moment of inertia about lateral body axis, slug- ft 2 (kg- m2) 
moment of inertia about normal body axis, slug-ft2 (kg- m2) 
product of inertia, slug-ft2 (kg-m2) 

radius of gyration in roll about principal longitudinal axis, ft (m) 

radius of gyration in yaw about principal normal axis, ft (m) 

nondimensional radius of gyration in roll about principal longitudinal axis, 
k Xo/ b 

nondimensional radius of gyration in yaw about principal normal axis, k^^ b 

vertical tail length, distance from moment reference center to aerodynamic 
center of vertical tail measured along fuselage center line, ft (m) 

airplane mass, slugs (kg) 

Mach number 

rolling moment, ft-lb (m-N) 
pitching moment, ft-lb (m-N) 
yawing moment, ft-lb (m-N) 
rolling velocity, rad/sec 
period of oscillation, sec 

effective dynamic pressure at vertical tail location, lb/ft2 (N/m2) 
free- stream dynamic pressure, lb/ f t^ (N/m2) 
yawing velocity, rad/sec 
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% 

t 


4/2 


wing area, ft 2 (m 2 ) 

nondimensional time parameter based on wing span, Vt/b 
time, sec 

time required for amplitude of oscillation to decrease by a factor of 

velocity, ft/sec (m/sec) 

angle of attack, deg or rad 

angle of sideslip, deg or .rad 

B .<b constants used in solution of characteristic equation 

^(r r O ,T 0 

X root of lateral-directional characteristic equation 

AX 4 + BX 3 + CX 2 + DX + E = 0 


V 


a 


0 


^b 

P 


<P 

* 


lateral-directional relative-density factor, m/pSb 
mass density of air, slugs/ft 2 (kg/m 3 ) 
sidewash angle, deg 
angle of bank, deg or rad 
angle of yaw, deg or rad 
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AC 


increment of C r 
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Model component designations: 


fuselage 


H 


horizontal tail 


vertical tail 


W 


wing 


ot 

due to vertical tail 


DESCRIPTION OF AIRPLANE 

The airplane studied in the investigation is a two-place twin- jet high-performance 
fighter designed for land and carrier-based operations. A three-view sketch showing the 
general layout of the configuration is presented in figure 1, geometric characteristics of 
the airplane are listed in table I, and typical mass characteristics for normal flight opera- 
tions (no external stores) are presented in table II. The longitudinal control system of the 
airplane consists of an all- movable horizontal tail (stabilator) which incorporates 23° 
negative dihedral (droop) to satisfy longitudinal stability requirements in the normal oper- 
ational flight range. The airplane lateral control system consists of upper-surface 
spoilers and ailerons. The control system is mechanized such that the ailerons deflect 
downward only while the spoilers deflect upward. The left aileron and right spoiler oper- 
ate simultaneously as do the right aileron and left spoiler. The directional control sys- 
tem consists of a conventional rudder. The maximum control- surface deflections are as 


follows: 

Rudder deflection, deg ±30 

Stabilator deflection (trailing edge), deg 21 up, 9 down 

Aileron deflection, deg 0 up, 30 down 

Spoiler deflection, deg 45 up, 0 down 


Mass loadings such as those presented in table II are typical of those of modern 
high-performance fighter airplanes and result from the fact that the major portion of the 
airplane mass is distributed along the fuselage; therefore the values of 1% and I Y are 
several times as great as those of lx- This type of inertial distribution has significant 
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effects on dynamic lateral-directional stability at high angles of attack as is discussed 
subsequently. 


FLIGHT MOTIONS 

The directional divergence exhibited by the airplane at high angles of attack is illus- 
trated by the time histories presented in figure 2. Shown in figure 2 are flight recorder 
traces of the major flight variables and control- surface deflections during an accelerated 
stall at 25 000 feet (7620 m) with the airplane configured for cruise flight (M = 0.4). 
Unfortunately there is no record of yaw angle or yaw rate. The maneuver was initiated 
by rolling to a 60° banked turn to the left. Angle of attack was then increased at an 
approximately constant rate. The normal acceleration trace indicates airframe buffet 
occurred at angles of attack as low as 10°. At about 38 seconds the magnitude of the 
normal acceleration trace starts to decrease, even though angle of attack is increasing, 
thereby indicating major stall. The angle of attack at this time was about 18°; this value 
should be remembered when analyzing the force test results presented subsequently. As 
the angle of attack increased further, lightly damped lateral oscillations about the longi- 
tudinal body axis (termed "wing rock") became noticeable. At about 44 seconds severe 
wing rock was experienced; at about 50 seconds the oscillation diverged violently and the 
airplane entered a 2^ -turn spin to the right. The flight path was about 40° below the 
horizon during the spin; therefore the spin appears as a continuous roll with reference to 
the earth axes. The angle of attack at the time of directional divergence was between 20° 
and 25°; this range of angle of attack should also be remembered for subsequent reference. 

METHOD OF ANALYSIS 

Static wind-tunnel force tests were conducted to determine the aerodynamic charac- 
teristics of the airplane at high angles of attack. Airframe components were tested indi- 
vidually and in several combinations to determine the contributions of the isolated com- 
ponents to the overall stability characteristics of the airplane and to determine mutual 
interference effects. Several geometric modifications or fixes to the basic configuration 
were evaluated in an attempt to delay or eliminate lateral-directional instability near the 
stall. 

The dynamic stability characteristics of the airplane were calculated by using linear 
three-degree-of-freedom equations of motion. The static and dynamic aerodynamic sta- 
bility derivatives used as input data in these calculations were measured quantities deter- 
mined in wind-tunnel tests of a model of the airplane. 
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Model and Test Equipment 

The wind-tunnel data presented herein were obtained with a 1/11 -scale model tested 
at the Langley Research Center in a low-speed wind tunnel with a 12-foot (3.66-m) octag- 
onal test section. Additional data obtained at higher values of Reynolds number during 
tests of a 1/15-scale model at the Ames Research Center are also presented for purposes 
of correlation. The two models differed in external engine inlet configuration and hori- 
zontal tail leading -edge configuration as shown in figure 3 - the inlets and horizontal tail 
of the 1/11-scale model being representative of those of an earlier airplane configuration 
and the inlets and horizontal tail of the 1/15-scale model being representative of those of 
a later airplane configuration. Both models had blocked inlets - that is, flow through the 
engines was not simulated. A number of geometric modifications to the basic airplane 
configuration were evaluated as possible fixes for the lateral-directional stability problem 
during the course of the study. These modifications, summarized in figure 4, consisted 
of vertical end plates on the stabilator, afterbody strakes, nose strakes, wing apex notches, 
wing fences, wing leading-edge droop, and two modified vertical tail surfaces. 

Tests 

Wind-tunnel force tests of the 1/11-scale model were conducted at low-subsonic 
speeds at a Reynolds number of 0.5 x 10® based on the mean aerodynamic chord of the 
wing. Measurements were made of the six force and moment components over an angle - 
of-attack range from 0° to 40° for a range of angle of sideslip of ±40°. The wing, fuselage, 
vertical tail, and horizontal tail were tested in several combinations. Tuft-flow visualiza- 
tion tests were conducted to aid in the interpretation and understanding of the force test 
results. A limited number of dynamic-pressure and sidewash measurements at the verti- 
cal tail location were also made in conjunction with an analysis of the effectiveness of the 
vertical tail at high angles of attack. Additional low-subsonic force test data measured 
at a Reynolds number of 4.3 x 10® at the Ames Research Center are also presented. 

Calculations 

The dynamic stability characteristics of the lateral -directional modes of motion 
were calculated for the basic and modified airplane configurations by using the linearized 
three-degree-of -freedom equations of motion presented in appendix A. The calculated 
characteristics included the period P and time to half- amplitude t ^ °* the Dutch 
roll mode and the time to half- amplitude of the spiral and roll subsidence modes. The 
magnitude and phase relationships of the roll <p, yaw if/, and sideslip /3 angles were 
also calculated to indicate the "shape" or general nature of the Dutch roll mode. Values 
of the dynamic lateral-directional stability derivatives used in the calculations were 
obtained from the results of unpublished forced-oscillation tests of the 1/11-scale model. 
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Dynamic stability calculations were made for trimmed level flight at an altitude of 
25 000 feet (7620 m). Mass characteristics used in the calculations are given in table II 
and the aerodynamic stability derivatives used are summarized in table m. 

RESULTS OF FORCE TESTS 

The present paper presents a general discussion of the results of all the pertinent 
force tests run during the study; however it presents only samples of the detailed data 
as needed to support the analysis. The results of the force tests are analyzed in terms 
of the individual contributions of the various airframe components to the overall aero- 
dynamic stability parameters. 

Longitudinal Stability Characteristics 

Basic configuration .- The static longitudinal characteristics of the basic configura- 
tion obtained during the force test program are presented in figure 5. The data denoted 
by a circular symbol are those obtained at low values of Reynolds number at the Langley 
Research Center; whereas the square symbols denote data obtained at considerably higher 
values of Reynolds number at the Ames Research Center. The data indicate good agree- 
ment between the high and low Reynolds number lift characteristics with the exception of 
the immediate poststall angle-of-attack range of 20° and 40°. These differences are 
believed to be due to the combined effects of Reynolds number and the geometric model 
differences mentioned previously. The variation of lift coefficient with angle of attack 
indicates that the lift-curve slope begins to decrease at an angle of attack of about 15°. 
Associated with the apparent flow separation is a reduction in the level of longitudinal 
stability. The low Reynolds number data indicate a condition of neutral longitudinal sta- 
bility from a = 15° to a = 20°; whereas the high Reynolds number data indicate a mild 
instability or pitch-up tendency for the same angle-of-attack range. This pitch-up ten- 
dency has been verified during stall entries with the full-scale airplane. The significance 
of the longitudinal aerodynamic characteristics in relation to the lateral-directional sta- 
bility of the airplane at the stall is twofold. First, the progression of major wing stall at 
an angle of attack of 15° should be reflected in changes in the lateral -directional aero- 
dynamic characteristics at the same angle of attack; and second, the existence of longi- 
tudinal instability at the stall would be expected to aggravate an already deteriorating 
control problem and cause the pilot to overshoot to the higher angles of attack associated 
with directional divergence. 

Modified configuration .- The effect of drooping the wing leading edge (40°) on longi- 
tudinal characteristics is shown in figure 6. It can be seen that the drooped leading edge 
maintained the linear portion of the lift curve up to a = 20° and increased the value of 
the maximum lift coefficient. The leading-edge droop also reduced the longitudinal 
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stability to such an extent that the configuration became neutrally stable over a much 
wider range of angle of attack. This reduction in stability is believed to be the result of 
leading-edge separation which is highly sensitive to variations in Reynolds number. 

Lateral-Directional Characteristics 

Basic configuration .- The variation of the static lateral-directional force and 
moment coefficients with angle of sideslip for the basic configuration are presented in 
figure 7. As can be seen, the variations of the lateral-directional coefficients were gen- 
erally linear over a sideslip range of ±20°. Stability derivatives obtained over a range 
of sideslip angle of ±5° are summarized in figure 8. The values of the stability deriva- 
tives are presented as the variations with angle of attack of the side-force derivative 
Cy^, the directional- stability derivative Cn^, and the effective-dihedral derivative C^. 
The circular symbols represent data obtained at a Reynolds number of 0.5 x 10®; whereas 
the square symbols denote data obtained at a Reynolds number of 4.3 x 10®. The two sets 
of data agree fairly well with the exception of a slightly lower level of directional stability 
for the higher Reynolds number data at moderate angles of attack. Both sets of data indi- 
cate a marked decrease in Cn^ as angle of attack is increased with the directional sta- 
bility being negative at angles of attack above 22°. The data also indicate that as the 
angle of attack exceeds 15° a substantial reduction in effective dihedral occurs, as might 
be expected from the lift curves presented in figure 5. Actually the loss of directional 
stability at a = 22° is considered to be especially serious because of the corresponding 
loss of effective dihedral. This aspect of the problem is discussed in detail subsequently. 
Also shown in figure 8 is a vertical hatched region indicating the maximum trim capabil- 
ity of the horizontal tail based on static force tests (i^ = -21°). It can be seen that there 
exists a trimmed angle-of-attack range of almost 10° for which the airplane is direction- 
ally unstable. 

Effect of airframe components .- The results of a series of tests to evaluate the 
contributions of various airframe components to the directional-stability derivative Cn^ 
are summarized in figure 9. Data obtained for the isolated wing and fuselage are pre- 
sented in figure 9(a) along with data obtained for the wing-fuselage combination. The 
data of figure 9(b) show the effects of the addition of the horizontal and vertical tail sur- 
faces to the wing-fuselage combination; and the data of figure 9(c) show the effects of the 
addition of the vertical tail to the isolated fuselage. Several points are to be noted from 
comparisons of the various data. For example, the isolated wing data W of figure 9(a) 
indicate little variation of Cn^ (relative to body axes) with increase in angle of attack. 
This result is to be expected inasmuch as little if any leading- edge suction is developed 
on thin sweptback wings. Also the isolated fuselage data F show that the fuselage does 
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not become more unstable until the angle of attack exceeds 28°. The wing-fuselage com- 
bination WF, however, shows large unstable variations of C n ^ with increasing angle of 
attack. The magnitudes of C n ^ are far greater than those expected by simple addition 
of the values for the isolated wing and fuselage. This type of aerodynamic phenomenon 
has been noted in the past for sweptback wing-fuselage combinations at high angle of 
attack (ref. 1). As stated in reference 1, the large unstable variations of directional 
stability which exist relative to the body axes are due largely to wing (or wing-fuselage) 
induced sidewash over the fuselage afterbody. An additional factor to be considered is 
the fact that when a swept-wing configuration is sideslipped at high angles of attack, 
leading-wing-panel stall may significantly affect the flow pattern at the rear of the air- 
plane. This fact is illustrated subsequently by tuft photographs. 

The adverse sidewash characteristics exhibited by the wing-fuselage combination 
might be expected to affect adversely the contributions of the vertical and horizontal tails 
to directional stability. The data of figure 9(b) support this assumption by showing that 
(1) the addition of the drooped horizontal tail H to the wing-fuselage combination WF pro- 
duced a large destabilizing effect at high angles of attack, and (2) the vertical tail con- 
tribution to directional stability diminishes markedly as a is increased above 20° and 
becomes destabilizing at angles of attack above 30°. 

The fact that the wing-fuselage combination, and not the fuselage itself, was respon- 
sible for the loss in tail effectiveness at high angles of attack is indicated by the data of 
figure 9(c) which show a comparison of the directional stability characteristics of the 
fuselage and vertical tail combination and the fuselage alone. This comparison shows 
that the vertical tail contribution at high angles of attack was reduced only slightly in 
comparison to the large loss in tail effectiveness shown by figure 9(b) when the tail was 
added to the wing-fuselage combination. The data of figure 9<c) show however that the 
effectiveness of the vertical tail at high angles of attack for the fuselage-tail combination 
was only about one-half as great as at low angles of attack. It is possible that this loss 
in tail effectiveness (note that the tail was never destabilizing) could have been partly or 
entirely the result of a reduction of dynamic pressure at the tail because of shielding by 
the fuselage at high angles of attack. 

In order to evaluate the adverse sidewash and reduction in dynamic pressure at the 
vertical tail location, measurements were made of the dynamic pressure and sidewash 
angle at the tail location. The results of the tests are shown as functions of angle of 
attack in figure 10. 

The results of the pressure survey are presented as qy/q^, the ratio of dynamic 
pressure at the vertical tail location divided by the free-stream dynamic pressure. The 
dynamic-pressure ratio is reduced as angle of attack is increased above 15°. More 
specifically, as the stall angle of attack is exceeded, the ratio is reduced to less than 
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50 percent of the value at a = 0° because of combined shielding of the vertical tail by 
the fuselage afterbody and impingement of the low velocity wake of the stalled wing- 
fuselage combination on the vertical tail. 

The sidewash angles at the vertical tail location were determined by comparing 
vertical tail-off data with data obtained by deflecting the vertical tail as an all- movable 
tail at various angles of sideslip relative to the fuselage. The results of the sidewash 

measurements are presented in terms of the sidewash parameter 1 - -§2 At low angles 

op 

of attack, the data of figure 10 indicate that the sidewash is quite favorable, as has been 
noted in the past for low -wing configurations. (See refs. 2 to 4.) At angles of attack near 
the stall, however, the sidewash parameter changes from favorable to unfavorable and, 
at a = 30°, the factor becomes negative, thereby indicating the existence of a destabil- 
izing flow field at the vertical tail location. 

The foregoing data and analysis show that the primary factor producing loss of 
directional stability at angles of attack immediately below the stall (angles in the 15° to 
20° range) is adverse sidewash on the fuselage afterbody and drooped horizontal tail 
induced by the wing-fuselage combination. As a is increased above a value of about 20°, 
the vertical tail enters the adverse sidewash field and the dynamic pressure at the verti- 
cal tail is reduced because of shielding by the aft fuselage and/or the stalled wake of the 
wing. Both of these factors contribute to a loss of vertical tail effectiveness at high 
angles of attack. The adverse sidewash effect is the predominant effect at angles of 
attack above 30° where the tail contribution to directional stability actually becomes 
negative. 

Modified configuration .- A number of geometric configuration modifications were 
tested in an effort to delay the loss of directional stability to higher angles of attack. 

Test results obtained with the modifications shown in figure 4(a) had no appreciable 
effect on either the angle of attack at which directional stability was lost or on the loss 
of effective dihedral above a = 15°. These results are therefore not presented in this 
paper. 

A second series of modifications was aimed at evaluation of the vertical tail size 
necessary to delay or eliminate the loss of directional stability. The results of tests 
conducted with the modified vertical tail surfaces shown in figure 4(c) are presented in 
figure 11. Increasing the aspect ratio or area of the basic vertical tail by 100 percent 
increased the directional stability at low angles of attack; however the configuration 
remained directionally unstable for angles of attack greater than 25°. It is interesting 
to note that increasing the vertical tail area made the configuration more unstable above 
a = 30° as would be expected based on the results of the sidewash measurements dis- 
cussed in the previous section. Increasing the aspect ratio did not produce this trend 
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probably because the increase in aspect ratio added vertical tail area above the adverse 
sidewash field. 

Inasmuch as it was determined that the wing-fuselage combination was a major 
factor leading to directional instability, tests were conducted to modify the basic-wing 
aerodynamic characteristics. It has been noted that the subject airplane appears to be 
more prone to exhibit directional divergence in the cruise configuration than in the 
landing-approach configuration (leading-edge slats, flaps, and gear extended). The air- 
plane in the landing- approach configuration can consistently be flown at the same angles 
of attack at which the cruise configuration exhibits violent lateral-directional instability. 
This difference is believed to be attributable to the effect of the extended slats and flaps 
on the aerodynamic characteristics of the wing-fuselage combination. Inasmuch as the 
full-scale airplane wing slats are equipped with boundary-layer control which is difficult 
to simulate at small scale, it was decided to obtain a gross indication of the effects of 
these devices by drooping the leading edge of the wing as shown in figure 4(b). The 
results of tests showing the effect of 40° leading-edge droop angle on Cn^ and C ^ are 
presented in figure 12. The data indicate that leading- edge droop produced two very sig- 
nificant and beneficial effects. First, the angle of attack at which directional stability 
was lost was increased from 22° to 25°; and second, was maintained at moderate 

levels up to a = 40°. These changes, especially the increase in at high angles of 

attack, have important beneficial effects on dynamic lateral-directional stability, as is 
shown subsequently herein. Additional tests were conducted to determine whether the 
increase in Cn^ at a = 25° was due to wing-fuselage characteristics or vertical tail 
effectiveness; the results of the tests are shown in figure 13. The increase in Cn^ 
appears to be related to both of these factors and is primarily an extension of vertical 
tail effectiveness to a = 25°. 

The previous results were obtained with both inner and outer panels of the wing 
leading edge deflected. (See fig. 4(b).) Several combinations involving various droop 
angles and leading-edge panels were tested; however the most effective combination was 
found to be 40° leading-edge droop on both inner and outer panels. 

RESULTS OF TUFT STUDIES 

A series of tuft studies were conducted to aid in the interpretation of the static force 
test results. Photographs of the tuft patterns as a is increased for sideslip angles of 
0° and -10° are presented in figure 14. At a = 0°, sideslip induces a strong favorable 
sidewash field at the vertical tail, as can be seen by the relative angles of the tufts on the 
fuselage afterbody. This result is in agreement with the results previously presented 
regarding variations of 1 - As the angle of attack is increased, however, the favor- 

able sidewash angle is reduced. At a = 25° and /3 = -10°, the tuft patterns indicate a 
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complete stall of the leading-wing panel and a reversed-flow region behind the stalled 
wing panel. This reversed-flow region is evident from the tufts on the fuselage behind 
the wing, and an extensive region of reversed or low-velocity flow in the vicinity of the 
vertical tail was found by exploring with a tuft on a wand. This reversed-flow and low- 
velocity region behind the stalled wing and fuselage is evidently the cause of the reduction 
of Qv/^oo at the rear °* the model, shown in figure 10, and is a major factor contributing 
to the directional instability of the wing-fuselage combination and to the loss in vertical 
tail effectiveness at high angles of attack. The results of the tuft studies also indicate 
that stall of the leading-edge panel is the cause of the loss of effective dihedral at 
a = 15°. Such a result is a characteristic of swept wings at high angles of attack. 

RESULTS OF DYNAMIC STABILITY CALCULATIONS 

The results of calculations to determine the lateral-directional dynamic stability 
of the basic configuration and of the modified configuration (incorporation of wing leading- 
edge droop) are presented in figure 15 and table HI. These results should be interpreted 
as qualitative (rather than quantitative) indications of the dynamic stability of the airplane 
for several reasons. First, the application of small perturbation, linearized equations of 
motion near the stall may be questionable because of the large amplitude motions involved 
(see fig. 2) and also because of the extremely nonlinear variations of yawing moment with 
angle of sideslip exhibited by the configuration at high angles of attack. (See fig. 7.) A 
second area of concern is that large variations in values of dynamic derivatives with 
amplitude and frequency at angles of attack near the stall have been noted in past inves- 
tigations. Also the calculations made for the modified configuration assumed changes in 
only the static stability derivatives Cn^, C^, and Cy^ and neglected possible varia- 
tions in dynamic derivatives because of wing leading-edge droop. The calculated results 
should however serve as a gross indication of the dynamic behavior of the vehicle for 
small disturbances from trimmed flight at large angles of attack. 

Basic Configuration 

The variation with angle of attack of the damping characteristics of the various 
modes of motion for the basic configuration are presented in the upper plot of figure 15(a). 
The damping of the roll, spiral, and Dutch roll modes are presented in terms of the 

damping factor —i— . Positive values of r-i— represent damped (dynamically stable) 
4/2 *1/2 

modes while negative values represent undamped (dynamically unstable) modes. The 
results of the calculations indicate that as a is increased, the roll and spiral modes of 
motion remain stable with the exception of the spiral mode at angles of attack between 18° 
and 22°. The unstable spiral mode indicated by the calculations has not however been 
evident in flight test results. The calculated results also indicate that the Dutch roll 
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oscillation becomes lightly damped at the stall and becomes unstable at an angle of attack 
of about 23°. This result agrees with the flight data presented in figure 2. The relative 
magnitudes and phase relationships of the roll-to-sideslip ratio <p/ j3 and roll-to-yaw 
ratio <p/\p presented in table III indicate the general nature of the lateral oscillation. 

The data indicate that as a is increased up to 15°, the magnitude of c p/ip is increased 
and the phase relationship of $/j3 approaches zero, thereby indicating that the Dutch 
roll becomes predominately a rolling motion about the longitudinal axis. This result is 
also in agreement with flight experience. At the poststall angles of attack (a = 25° 
and 30°), the phase and magnitude relationships change such that the magnitude of <p/xp 
is considerably reduced and 6 and \p become approximately in phase. These char- 
acteristics are indicative of larger yawing content in the Dutch roll mode which would, 
because of the relative phasing of <p and \p and the high degree of instability (time 
required to double amplitude less than 2 seconds), appear to the pilot as a yawing or 
directional divergence. These results indicate that the directional divergence exhibited 
by the airplane is in reality a highly unstable Dutch roll oscillation, the character of which 
changes from that of a predominately rolling motion at prestall angles of attack to that of 
a predominately yawing motion at poststall angles of attack. 

In the past, numerous investigations of lateral-directional dynamic stability at high 
angles of attack have been conducted. (For example, see refs. 5 to 9.) One common and 
significant outcome of the studies was the indication that configurations having relatively 
high values of the ratio I z /lx together with appreciable values of the effective-dihedral 
parameter -C^ usually diverge under dynamic flight conditions at angles of attack 
higher than those expected based on the variation of Cn^ with a. This characteristic 
is attributed to the favorable effect of -C ^ on dynamic stability as described in appen- 
dix B. As a result of past experience the parameter C n ^ dynamic is usually used as 

an indication of directional divergence at high angles of attack. Shown in the lower plot 
of figure 15(a) are the variations of C n ^ and C n ^ dynamic with angle of attack for the 

basic configuration. The data show that C n ^ and C n ^ dynamic a PP roac ^ zero at a 

common angle of attack of about 21°. This result is attributable to the reduction of -C ^ 
at the stall as noted in figure 8. The angle of attack at which Cn^ dynamic t> ecomes 

negative does give a good approximation of the angle of attack at which the damping of the 
Dutch roll mode becomes unstable. 

Modified Configuration 

The results of the calculations for the airplane with 40° leading-edge droop are 
presented only to illustrate the effects of leading-edge devices on lateral-directional 
stability. Leading-edge droop is not proposed herein as an acceptable modification to the 
airplane because of the reduction of longitudinal stability shown in figure 6. 
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The variation of the damping of the lateral-directional modes of motion for the 
modified configuration are shown in the upper plot of figure 15(b). The data indicate that 
all modes of motion remain dynamically stable up to a = 30°. The elimination of 
dynamic instability is produced by the increases in C n ^ and given in figure 12. 


The lower plot in figure 15(b) indicates that C 


n /3, dynamic 


also remained positive up to 


a = 30° because of the moderate level of C ^ afforded by leading- edge droop. 


The results obtained with leading-edge droop indicate that leading-edge devices can 
significantly alter the dynamic lateral-directional behavior of the present configuration. 
The effects of leading-edge droop are believed to be similar to those obtained with the 
airplane configured for the power approach. The foregoing analysis explains why the 
airplane remains dynamically stable in the approach configuration. It should also be 
emphasized that the gains afforded by high-lift devices are limited and that the airplane 
would become unstable at higher angles of attack. 


CONCLUSIONS 


The results of an analysis of the lateral-directional stability characteristics of a 
twin-jet swept-wing fighter airplane at high angles of attack indicate the following 
conclusions: 

1. The directional divergence exhibited by the airplane is brought about by a simul- 
taneous loss of directional stability and effective dihedral at high angles of attack. 

2. The loss of directional stability results from a combination of an adverse side- 
wash region at the rear of the airplane and a reduced dynamic pressure at the vertical 
tail location. The adverse sidewash is generated by the wing-fuselage combination and 
is related to stalling of the leading wing during sideslip at high angles of attack. 

3. The loss of effective dihedral is a result of stalling of the leading-wing panel at 
high angles of attack. 

4. The apparent directional divergence is in reality part of one cycle of a highly 
unstable Dutch roll oscillation. 

5. The only geometric modification studied that significantly delayed the divergence 
was wing leading- edge droop. The leading- edge droop however reduced the level of lon- 
gitudinal stability for the low value of Reynolds number used in the tests. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., June 4, 1969, 

126-62-01-01-23. 
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APPENDIX A 


EQUATIONS OF MOTION 

The nondimensional lateral-directional equations of motion referred to a principal 
body-axis system are as follows: 

Roll: 

- 2 Mb K X o 2D b 2 0 + + \ C ip D b^ + 2 ~ 0 

Yaw: 

-2fX b K Zo 2 D b V + Cnp/3 + ± C np D b 4> + i = 0 

Sideslip: 

-2M b (D b /3 + cos + C Y/3 0 + \ c Y p D b 4> + 2 M b sin o;D b <^> 

+ Cl sin a\}/ + Cl cos a<p + ^ = 0 

\Si~ XSi^ XSr^ 

When <p Q e D is substituted for <p, \J/ Q e for \p, and /3 Q e for 0 in the 
equations written in determinant form, X must be a root of the characteristic equation 
(neglecting the zero or heading root) 

AX 4 + BX 3 + CX 2 + DX + E = 0 


where 


A = 8|i b 3 K Xo 2 K Z() 2 

B = -2n b 2|2K Xo 2 K Zo 2 C Yj3 * K Xo 2 C„ r + K Z() 2 C ;p ) 

C * “b(- K Xo 2c Y r S + % 0 2c ”r C Y/3 ♦ ^b’V’S C0S “ + K Z 0 2c ip% 
- K Zo 2 %C Yp - sin or + | C„ r C ip - i C„ p C^ 
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APPENDIX A 


D = -M b ( 2K Xo 2 C n ^C L sin or + 2K Zo 2 C^C L cos a + C n ^C lr sin a - C^Cn r sin a 
+ C n pC lp cos a - C Z(3 C np cos aj - ^Cn r Cj p C Y 0 + Cn r C i(3 C Yp - c n p c Y r c £/3 
+ CnpC^ r C Y ^ - C n ^C Yp q r + Cn^C Yr C^ p ^ 

E = |c L cos a(c nr C lfi - C lr C^ + | C L sin <*(c lp C nj3 - C np C^ 


The period and damping of a mode of motion in seconds are given, respectively, by 
the equations 

p = 2 jt b 
d V 


and 

* _ 0.693 b 

*1/2 ~ c y 


where c and d are the real and imaginary parts of the root of the characteristic 
equation. 



APPENDIX B 


INTERPRETATION OF <fy fdynftmlc 

During the course of the study, it was observed that the lateral-directional diver- 
gence occurring in the poststall region was not an aperiodic divergence but rather a 
violently unstable oscillation. Inasmuch as the parameter 

I z 

dynamic = Cn |3 ' % sin “ 

has been shown to correlate with directional divergence of inertially slender configura- 
tions, an attempt was made to correlate the calculated oscillatory instability with the 
parameter C n ^ dynamic' K was observed that the oscillatory instability was associated 

with a change in sign of the coefficient C of the characteristic equation. (See appen- 
dix A.) This particular coefficient contains several combinations of derivatives; how- 
ever most of the magnitude of the coefficient is composed of the term 

4 M b (K X() 2 C n p cos a - K Zq 2 C ^ sin a) 

This expression C n ^ dynamic > ® therefore is an approximate criterion for diver- 
gence in the form of lateral-directional oscillatory instability. 
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TABLE I.- DIMENSIONAL CHARACTERISTICS OF AIRPLANE 


Overall length 

Wing: 

Span 

Area (including leading-edge extension) . . . 

Root chord 

Tip chord 

Mean aerodynamic chord, c 

Leading edge of c rearward of leading edge 

of root chord 

Aspect ratio 

Taper ratio 

Sweepback of 25 percent chord 

Dihedral (inboard 69.5 percent b/2) 

Dihedral (outboard 69.5 percent b/2) . . . . 

Incidence 

Airfoil section: 

Root 

Tip 

Aileron: 

Area (one side) rearward of hinge line . . 
Span (one aileron) (from 44.5 percent b/2 


to 67.0 percent b/2) 4.35 ft 

Inboard end chord (base line 

103.24 in. (262.23 cm)) 37.81 in. 

Outboard end chord (base line 

155.44 in. (394.82 cm)) 34.38 in. 

Spoilers: 

Area (one side) 


Span (from 45.3 percent b/2 to 67.0 percent b/2) 

Inboard end chord 

Outboard end chord 

Horizontal tail: 

Area (in chord plane) 

Movable area 

Span 

Aspect ratio 

Taper ratio 

Sweepback of 25-percent chord 


57.59 ft (17.55 m) 

.... 38.41 ft (11.71 m) 

. . . 538.34 ft 2 (50.01 m 2 ) 
. . 282.00 in. (716.28 cm) 
. . 47.00 in. (119.38 cm) 

. . 192.50 in. (488.95 cm) 

. . 110.76 in. (281.33 cm) 

2.82 

. 0.167 

45.00° 

0 ° 

12 . 00 ° 

1 . 00 ° 

NACA 0006.4-64 (modified) 
NACA 0003.0-64 (modified) 

13.08 ft 2 (1.22 m 2 ) 

(1.33 m) (22.5 percent b/2) 


(96.04 cm) (21.3 percent c) 
(87.33 cm) (27.6 percent c) 


13.08 ft 2 (1.22 m 2 ) 
4.19 ft (1.28 m) 
1.39 ft (0.42 m) 
1.04 ft (0.32 m) 

94.70 ft 2 (8.80 m2) 
77.40 ft 2 (7.19 m2) 
17.705 ft (5.40 m) 

3.30 

0.20 

35.50° 
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TABLE I.- DIMENSIONAL CHARACTERISTICS OF AIRPLANE - Concluded 


Dihedral 

Root chord (at airplane center line) 

Tip chord (theoretical) 

Airfoil section: 

Root (airplane center line) 

Tip (theoretical) 

Hinge- line location, percent c^ 

Vertical tail: 

Area 

Span 

Taper ratio 

Root chord , 

Tip chord ...... 

Sweepback of 2 5 -percent chord 

Airfoil section: 

Root 

Tip 

Rudder: 

Area (rearward of hinge line) 

Hinge-line location, percent of water-line chords 


-23.00° 

107.00 in. (271.78 cm) 
21.40 in. (54.36 cm) 


NACA 0003.7-64 (modified) 
NACA 0003.0-64 (modified) 
41.00 

. . . . 67.50 ft 2 (6.27 m2) 

. . . . 6.38 ft (1.94 m) 

0.227 

. . 207.15 in. (526.16 cm) 
. . 47.10 in. (119.63 cm) 

58.30° 

NACA 0004.0-64 (modified) 
NACA 0002.5-64 (modified) 

. . . . 11.07 ft 2 (1.03 m 2 ) 
80.00 
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TABLE H.- MASS AND INERTIA CHARACTERISTICS OF AIRPLANE 


Weight, lb (N) . . . . 39 099 (173 912) 

Center-of-gravity position (percent c) 33 

Relative density, jji b : 

Sea level 24.73 

25 000 ft (7620 m) 55.16 

Moments of inertia: 

I x , slug-ft 2 (kg-m 2 ) 29 950 (35 183) 

I Z , slug-ft 2 (kg-m 2 ) 169 538 (229 860) 

IXZ’ slu g- ft2 (kg-m 2 ) 5241 (7106) 
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TABLE III.- CALCULATED DYNAMIC STABILITY CHARACTERISTICS 
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'Spiral and roll modes combine to form second lateral oscillation having t^ = 0-91 sec and P = 18.7 sec. 
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Figure 2.- Concluded. 
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Figure 3.- Sketches of geometric differences between models. 










Figure 4.- Continued. 
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Figure 5.- Variation of static longitudinal characteristics with angle of attack, if = 0°. 
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Figure 6.- Effect of wing leading-edge droop on longitudinal characteristics, ij = 0°. 


32 



1.0 



Figure 7.- Variation of static lateral-directional characteristics with angle of sideslip. 
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Figure 7.- Concluded. 


34 



a, deg 


Figure 8.- Variation of static lateral-directional stability derivatives with angle of attack. 


35 









O Basic vertical tail 
A Aspect ratio increased 100% 
O Area increased 100% 

□ Vertical tail off 



Figure 11.* Effect of vertical tail modifications on C nR - 
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(a) Basic configuration. 

Figure 15.* Results of dynamic stability calculations. 
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